Bound states in continuum in an electron waveguide 
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It is shown that in a double-cavity, two-dimensional electron waveguide, the interplay between 
quasi-bound states of each cavity leads to the appearance of bound states in continuum for certain 
distances between the cavities. These bound states may be used to trap electrons in de-localized 
states distributed in both cavities. 
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Ever since von Neumann and Wigner 0] proposed that 
certain type of oscillating attractive potentials could pro- 
duce isolated bound states with energies within the con- 
tinuum 0, El- a number of studies have been reported 
the presence of the "bound states in continuum" (BIC) 
that can exist above the continuum minimum. Fonda 
and Newton discussed BIC in a system of two coupled 
square-well potentials using resonance scattering theory 
Trajectories of the poles of this system have been 
studied in Refs. Positive energy bound states in 

super-lattice structures with a single impurity potential 
or a single defect stub @ have been reported. 

In Ref. [9| the existence of BIC was demonstrated 
for a model of two atoms in one-dimension. The BIC 
were formed by the exchange of resonant photons. A 
physically analogous effect occurs in a two-dimensional, 
double cavity waveguide shown in Fig. ^ The purpose of 
the present study is to show the the existence of BIC in 
this waveguide configuration. Electron waveguides can 
be formed at a GaAs/AlGaAs interface [Icj . 

The leads of the waveguide in Fig. ^ form an infinite 
quasi-one-dimensional wire, where electrons have a con- 
tinuous spectrum of energy. Due to the lateral confine- 
ment in the wire, the spectrum has a minimum energy 
that allows propagation along the wire. If there is a single 
cavity, an electron inside the cavity with energy below the 
minimum will form bound states. In contrast, an electron 
with energy above the minimum will form quasi-bound 
states with finite life-time, where the electron escapes the 
cavity through the leads. Bound states and quasi-bound 
states are associated with real and complex energy poles 
of the scattering matrix . 

For the double-cavity waveguide, the electron states 
are much more varied. In this study, we focus on the 
interplay between the quasi-bound states formed in each 
cavity. We find BIC appear for certain distances between 
the cavities. These positive energy bound states are quite 
different from previous ones since they are formed by the 
interaction between two quasi-bound states. 

Theoretical model. We use the single-electron Hamil- 
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FIG. 1: Double cavity electron waveguide 
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with vanishing wavefunctions at the boundaries. The ef- 
fective mass m* is m* = 0.05to c , where m e is the mass of 
the free electron. We assume the two cavities are iden- 
tical. To analyze this Hamiltonian, we decompose the 
wave guid e into two independent closed cavities and the 
lead • The wavefunctions inside the cavities are de- 
noted by 1 771, n)i, where i = 1,2 labels the cavities and 
to, n are positive integers representing the horizontal and 
vertical wave numbers. The wavefunctions in the leads 
are denoted by |fc, j) where k is the horizontal wave num- 
ber (real) and j the vertical wave number (integer). The 
energies of the wavefunctions in either cavity and the lead 
are, respectively 
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The minimum energy for propagation along the lead is 
Ei(0, 1). We consider an electron with low energy nar- 
rowly centered around 



E° = E c (m ,n ) 



(3) 



We assume that Ei(0, 1) < E° c < Ei(0,2). The electron 
may propagate through the first mode of the lead, but 
not through the higher (j > 1) modes. In our theoretical 
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model, we neglect the j > 1 modes, which are evanescent, 
and keep only the j = 1 mode. Henceforth we omit the 
j = 1 index, e.g., Ei(k) = Ei(k, 1). 

We will rewrite the Hamiltonian using the cavity and 
lead basis states. With \i) — \m ,no)i, we define the 
cavity and lead projectors 
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To make the basis states orthogonal, we introduce the 
modified lead states 

|^ fc ) = \k) -RP c (P c RP c )~ 1 \k), 

R = Pl E l{ k)-H + l Pl (5) 

that satisfy (i\ip k ) = and (Vv \H \ip k ) = Ei(k)S(k - k'). 
The following approximate Hamiltonian is obtained 
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The terms Vi(k) — (i\H\ipk) represent the amplitude of a 
transition of the electron from the lead to the cavities or 
vice versa. For cavities centered at x = x\ and x — X2, 
they have the form 



V h2 (k)^v k e lkx ^ +u k e l 
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To construct eigenstates of the Hamiltonian JBJ we start 
with the symmetric and anti-symmetric states 



|±) = (|1)±|2))/V2 



(8) 



From these we obtain eigenstates with complex energy 
eigenvalues z° Q , which are poles of the S-matrix. For 
t > we take the eigenvalues with negative imaginary 
part. They are solutions the integral equation 
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closest to the real axis. Here d = \x% — xt\ is the distance 
between the cavities. The + superscript means analytic 
continuation from the upper to the lower half-plane of 
z± . As the distance d is varied, the poles z± move in the 
complex plane. For certain distances d± the imaginary 
part of z± vanishes || . This happens when 1 ± cos kd± = 
for Ei(k) = z±. These conditions give 



2n+(l±l)/2 
d± = — -. Trn 



(10) 



with n integer. Replacing d = d± in Eq. @ we obtain 
real solutions for z±, which suggests the existence of BIC. 
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FIG. 2: Transmission probability, |T| 2 , versus energy, E, for 
the double-cavity waveguide with W c = 2, L c = 2, Wi = 1, 
and (a) d = 5.25, (b) d = 5.5, (c) d = 5.60, (d) d = 5.70, 
and (e) d = 5.90. The dashed line in (c) is the transmission 
probability for the single-cavity waveguide. 



In the following we verify this using a more accurate de- 
scription of the electron waveguide. 

Computational results. We have computed the exact 
energy eigenstates for the lowest propagating mode in 
the double-cavity waveguide as a function of energy and 
the distance between the two cavities using the bound- 
ary integral method [ill fl3| . The exact eigenstates are 
built out of local propagating and evanescent modes in 
the leads and cavities and are composed of incoming, 
ijj~{x,y), and outgoing, ip + (x,y) states, 



i>{n,y) = i> (x,y) + S(E)ip + (x,y), 



(11) 



where the scattering amplitude, S(E), is composed of 
reflection and transmission coefficients. The unit of the 
length is the width of the lead, Wi = 1, which corre- 
sponds to lOOA. The width and length of each cavity is 
W c = 2 and L c — 2, respectively. 

We have computed the transmission probability, |T| 2 , 
for the lowest propagating mode as a function of energy 
for the double-cavity waveguide in Fig. [21 Transmission 
zeros in a given channel have been related to Fano type 
resonances between the continuum of that channel and 
bound states of the closed cavity of higher energy 0, 
1 1 5| . The strength of the coupling between the continuum 
state and the bound state determines the decay rate and 
the width of the transmission profile 0, 0] . 

As we vary the distance between the two cavities, the 
transmission profiles change near the resonance energy 
region E — 0.25eV. There are sharp peaks in the trans- 
mission profiles on either side of the resonance energy 
regions in Figs. [3(b) and (d). There is a broad transmis- 
sion profile with d = 5.60 in Fig. [2(c). For comparison, 
the transmission probability for the single-cavity waveg- 
uide is shown in Fig. El(c) as a dotted line. The transmis- 
sion zero profile shows a wider dip for the double-cavity 
waveguide compared with that of the single-cavity waveg- 
uide. We will see later that the BIC is formed when the 
distance between the two cavities is d = 5.60. These fea- 
tures are associated with the presence of the double poles 
in the complex energy plane and affect the dynamics of 
the electron in the double-cavity waveguide. 

The transmission zeros in Fig. [21 are associated with 
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FIG. 3: Transmission amplitude, T(E), in the complex energy 
plane, -E = Er + iEi, for the double- cavity waveguide with 
W c = 2, L c = 2, W; = 1, and (a) d = 5.25, (b) d = 5.5, (c) 
d = 5.60, (d) d = 5.70, and (e) d = 5.90. The transmission 
amplitude in the complex energy plane for the-single cavity 
waveguide is in (f). The dark regions indicate the positions of 
the poles and the bright region indicates the position of the 
transmission zero. 
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FIG. 4: Distribution of the pole pairs in the complex energy 
plane as we vary the distance between the two cavities. The 
dark point on the real axis indicates the energy of the trans- 
mission zero. 



the poles of T(E) in the complex energy plane. Transmis- 
sion amplitude in the complex energy plane has a branch 
cut starting from the lower edge of the continuum and 
extending along the positive energy axis and has poles 
at energy values Er — i-y. These poles give rise to the 
transmission zeros on the positive real axis [TlL Il7| . 

The locations of the poles in the complex energy plane 
for the double-cavity waveguide are shown in Figure 03 
As we increase the distance between the two cavities, 
the pole located on the right-hand side in Fig. 03 (a) 
approaches the real axis. The pole disappears into the 
real axis and produces a zero value of 7 for the distance, 
d = 5.60, in Fig. 03 (c). This implies the formation of 
a BIC with infinite life-time. The pole on the left-hand 
side in Fig. 03 (a) moves away from the real axis and 



attains large decay rate 7, accordingly. As the distance 
between the cavities is further increased, the pole on the 
real axis in Fig. 03 (c) emerges out of the transmission 
zero and recedes from the real axis. If we continue to 
increase d, the pole on the right-hand side in Fig. 03 
(e) approaches the real axis and makes a transit through 
the transmission zero at d = 6.26 (not shown). BIC 
states regularly reappear as the distance between the two 
cavities is varied. This agrees with the existence of real 
solutions of Eq. @ with Eq. I|10|) for different values of n. 
We have found BIC states for distances of up to d ~ 50, 
which means that these states can be quite de-localized. 

In a separate calculation, the wavefunctions inside the 
two cavities for the real energy states associated with the 
right-hand side and left-hand side poles in Fig. 03(a) show 
symmetric and anti-symmetric structures, respectively. 
These wavefunctions are similar to the symmetric and 
anti-symmetric combinations of the eigenstates of the two 
closed cavities, 

tp±(x,y) = -i= (sin(27r(x- xi)/i c )cos(37ry/W / c )xi 

± sm(2ir(x-x 2 )/L c )cos(3ny/W c )x2) (12) 

where Xi — \ inside cavity i and \i = outside. Eq. i|12|) 
corresponds to Eq. ijSJl • 

Each complex pole keeps its symmetric or anti- 
symmetric feature as it migrates throughout the com- 
plex energy plane, as we vary the distance between the 
two cavities. Thus the BIC sstates of the electron ap- 
pearing at d = 5.6 and d — 6.26 are symmetric and anti- 
symmetric states, respectively. Solutions of Eq. © for 
n = A and n = 5 give d+ = 6 and d_ = 6.67, respectively, 
which are in fair agreement with the computational re- 
sults. 

In Figure 01 we show the distribution of the pole pairs 
in the complex energy plane as we vary the distance be- 
tween the two cavities. Each of the poles with either 
symmetric or anti-symmetric identity makes a counter- 
clockwise circulation that passes through the transmis- 
sion zero on the real axis. The dark point on the real 
axis indicates the energy of the transmission zero. It has 
been known that the double or multiple poles induced 
by a laser in an atom coalesce to form an exceptional 
point or repel each other to form an avoided crossing 
|l8l flii. Eof . In contrast to these works, the poles of the 
double-cavity waveguide make a circular motion and do 
not approach or repel each other. 

One of the important features of an electron waveg- 
uide is that the energy spectrum of the electrons is con- 
tinuous. This means that the electron waveguides fall in 
the class of unstable systems. An electron in the cavity 
with wavelength smaller than the critical wavelength for 
propagation along the lead, will ultimately decay (escape 
through the leads). However, in our study, we have shown 
that some of the states of an electron in a double-cavity 
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FIG. 5: The survival probability, P(t) — \A(t)\ 2 , versus time, 
t for the states prepared symmetrically or anti-symmetrically 
inside the two cavities at t = 0. The distances between the 
two cavities are (a) d — 5.60 and (b) d = 6.26. The dotted line 
shows the survival probability of a state in the single-cavity 
waveguide in (a) and (b). The solid line displays the survival 
probability of a state prepared symmetrically and the dashed 
line displays the survival probability of a state prepared anti- 
symmetrically. 



In conclusion, we have proved that there can be BIC 
states with specially arranged geometry in double-cavity 
electron waveguides. These maybe used as a quantum in- 
formation storing device. Pairs of electrons with opposite 
spins can form de-localized, entangled states. The exis- 
tence of BIC states may be verified experimentally using 
actual electron waveguides. Alternative experimental se- 
tups are electromagnetic waveguides, which are described 
by a similar model, and semiconductor arrays of quantum 
wells H3. 
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waveguide can prevent the state from decaying and trap 
the electron inside the two cavities with infinite life-time, 
forming BIC. 

The survival probabilities of an electron placed in the 
waveguide cavities can be calculated using the scatter- 
ing states of the double-cavity waveguide. The survival 
probability can be written as P^(t) — \A x p{t)\ 2 , where 
A^(t) is the survival amplitude, 

poo 

A^(t) = (V>|e-^*/ s |V> = / dE\(i,\E)\ 2 e~ lEt ' h . (13) 

Jo 

We choose as initial state the symmetric or anti- 
symmetric combinations of the eigenstates of the two 
closed cavities in Eq. (|12[) . with no probability amplitude 
outside the cavities. In order to solve Eq. (|13|) numeri- 
cally, we discretize the energy eigenstates, \E), residing 
in the continuum. 

In Fig. EJ we plot the survival probability, P{t) = 
\A(t)\ 2 , versus time, t for the states prepared symmet- 
rically or anti-symmetrically inside the two cavities at 
t = 0. The distances between the two cavities are 
d = 5.60 and d — 6.26 for Figs. El (a) and (b), respec- 
tively. With this arrangement, one of the complex poles 
has a disappearing imaginary part, 7 — + 0, and gives rise 
to an infinite life-time. The results show that the sym- 
metrically (anti-symmetrically) arranged state does not 
decay over time for the cases with d = 5.60 (6.26). On 
the other hand, the anti-symmetric (symmetric) state de- 
cays quickly. The dotted line in Fig. Elshows the survival 
probability of the state, ip — sin(27ra; / L c )cos(3ny / D c ) , in 
the single-cavity waveguide. 
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